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We examine the efficacy of streamwise traveling waves generated by a zero-net-mass-flux 
surface blowing and suction for controlling the onset of turbulence in a channel flow. 
For small amplitude actuation, we utilize weakly nonlinear analysis to determine base 
flow modifications and to assess the resulting net power balance. Receptivity analysis of 
the velocity fluctuations around this base flow is then employed to design the traveling 
waves. Our simulation-free approach reveals that, relative to the flow with no control, 
the downstream traveling waves with properly designed speed and frequency can sig- 
nificantly reduce receptivity which makes them well-suited for controlling the onset of 
turbulence. In contrast, the velocity fluctuations around the upstream traveling waves 
exhibit larger receptivity to disturbances. Our theoretical predictions, obtained by per- 
turbation analysis (in the wave amplitude) of the linearized Navier-Stokes equations with 
spatially periodic coefficients, are verified using full-scale simulations of the nonlinear flow 



dynamics in companion paper, Lieu, Moarref & Jovanovic (2010). 



1. Introduction 

The problem of turbulence suppression in a channel flow using feedback control with 
wall-mounted arrays of sensors and actuators has recently received a significant atten- 
tion. This problem is viewed as a benchmark for turbulence suppression in a variety of 
geometries, including boundary layers. Also, there has been mounting evidence that the 
linearized Navier-Stokes (NS) equations represent a good control-oriented model for the 
dynamics of transition. Recent research suggests that, in wall-bounded shear flows, one 
must account for modeling imperfections in the linearized NS equations since they are 
exceedingly sensitive to external excitations and unmodelled dynamics (see, for example, 
Trefethen efar|l993| |Farrell fc Ioannou||1993| |Jovanovic fc Bamieh|2005| |Schmid|2007[ ) . 



This has motivated several research groups to use the linearized NS equations for model- 



based design of estimators and controllers in a channel flow (Bewley & Liu 1998; Lee 
et a/.|2001|rKim|2003|[Hogberg et a/.|2003a|fr||Hcepffner et a/.|2005||Chevalier et a/.|2006 



Kim fc Bewley||2007|)Vazquez fc Krstic||2007a|6| |Cochran fc Krstic||2QQ9| ). These results 



suggest that the proper turbulence suppression design paradigm is that of disturbance 
attenuation or robust stabilization rather than modal stabilization. 

An alternative approach to feedback flow control relies on the understanding of the ba- 
sic flow physics and the open- loop implementation of controls (i.e., without measurement 
of the relevant flow quantities and disturbances). Examples of sensorless strategies in- 
clude: wall geometry deformation such as riblets, transverse wall oscillations, and control 
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of conductive fluids using the Lorentz force. Although several numerical and experimen- 
tal studies show that properly designed sensorless strategies may yield significant drag 
reduction, an obstacle to fully utilizing these physics-based approaches is the absence of 
a theoretical framework for their design and optimization. 

An enormous potential of sensorless strategies was exemplified by |Min et al. (2006), 
where direct numerical simulations (DNS) were used to show that a surface blowing and 
suction in the form of an upstream traveling wave (UTW) results in a sustained sub- 
laminar drag in a fully developed turbulent channel flow. The underlying mechanism 
for obtaining drag smaller than in a laminar flow is the generation of the wall region 
Reynolds shear stresses of the opposite signs compared to what is expected based on 
the mean shear. By assuming that a wall actuation only influences the velocity fluc- 
tuations, Min et al. (2006) determined an explicit solution to the two dimensional NS 
equations linearized around parabolic profile; they further used an expression for skin- 



friction drag in fully developed channel flows (Fukagata, Iwamoto & Kasagi 2002 Bewley 



fc Aamo||2004| ), and showed that the drag is increased with the downstream traveling 



waves (DTWs) and decreased with the upstream traveling waves. 

A comparison of laminar and turbulent channel flows with and without control was 
presented by Marusic, Joseph & Mahesh (2007), where a criterion for achieving sub- 



laminar drag was derived. This study considered effectiveness of streamwise traveling 
waves at high Reynolds numbers and discussed why such controls can achieve sub-laminar 
drag. Another recent study, Hcepffner fc Fukagata| ( |2009[ ) , emphasized that the UTWs 
introduce a larger flux compared to the uncontrolled flow which motivated the authors 
to characterize the observed mechanism as a pumping rather than as a drag reduction. It 
was shown that, even with no driving pressure gradient, blowing and suction along the 
walls induces pumping action in a direction opposite to that of the wave propagation. 
By considering flows in the absence of velocity fluctuations Hcepffner fc Fukagata| (2009) 
showed that it costs more to drive a fixed flux with wall-transpiration type of actuation 
than with standard pressure gradient type of actuation. A fundamental limitation on 
the balance of power in a channel flow was recently examined by |Bewley| ( |2009[ ); this 
study showed that any transpiration-based control strategy that results in a sub-laminar 
drag necessarily has negative net efficiency compared to the laminar flow with no control. 
Furthermore, Fukagata, Sugiyama & Kasagi (2009j) showed that a lower bound on the net 
driving power in a duct flow with arbitrary constant streamline curvature is determined 
by the power required to drive the Stokes flow. It was thus concluded that the flow has 
to be relaminarized in order to be driven with the smallest net power. However, since 
the difference between the turbulent and laminar drag coefficients grows quadratically 
with the Reynolds number, Marusic et al. (2007) argued that relaminarization may not 
be possible in strongly inertial flows. An alternative approach is to design a controller 
that reduces skin-friction drag in turbulent flows; provided that the control power is less 
than the saved power, a positive net efficiency can still be achieved. 

In this paper, we show that a positive net efficiency can be achieved in a channel 
flow subject to streamwise traveling waves if the controlled flow stays laminar while 
the uncontrolled flow becomes turbulent. Starting from this observation, we develop a 
framework for design of the traveling waves that are capable of (i) improving dynamical 
properties of the flow; and (ii) achieving positive net efficiency. We quantify receptivity 
of the NS equations linearized around UTWs and DTWs to stochastic disturbances by 
computing the ensemble average energy density of the statistical steady-state. Motivated 
by our desire to have low cost of control we confine our study to small amplitude blow- 
ing and suction along the walls. This also facilitates derivation of an explicit formula 
for energy amplification (in flows with control) using perturbation analysis techniques. 
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Our simulation-free design reveals that the UTWs are poor candidates for preventing 
transition; conversely, we demonstrate that properly designed DTWs are capable of sub- 
stantially reducing receptivity of three dimensional fluctuations (including streamwise 
streaks and Tollmien-Schlichting (TS) waves). This indicates that the DTWs can be 
used as an effective means for controlling the onset of turbulence. Moreover, we show the 
existence of DTWs that result in a positive net efficiency compared to the uncontrolled 
flow that becomes turbulent. Our theoretical predictions are verified in Part 2 of this 
paper ( [Lieu et aZ.|2QlQ| ) using DNS of the NS equations. Thus, our work (i) demonstrates 



that the theory developed for the linearized equations with uncertainty has considerable 
ability to capture full-scale phenomena; and (ii) exhibits the predictive power of the 
proposed perturbation- analysis-based method for designing traveling waves. 

This paper represents an outgrowth of the study performed during the 2006 Center 
for Turbulence Research Summer Program ( |Jovanovic, Moarref fc Y ou 2006). While [ja 



vanovic et al\ (|2006j) only focused on receptivity of UTWs with large wavelength, our 
current study does a comprehensive analysis of the influence of both UTWs and DTWs 
on the fluctuations' kinetic energy and the overall efficiency. We also note that linear 



stability and transient growth of traveling waves were recently examined by |Lee, Min 



& Kim (2008). For selected values of parameters, it was shown that the UTWs destabi- 
lize the laminar flow for control amplitudes as small as 1.5% of the centerline velocity; 
on the other hand, the DTWs with phase speeds larger than the centerline velocity re- 
main stable even for large wave amplitudes. Moreover, the UTWs (DTWs) exhibit larger 
(smaller) transient growth relative to the uncontrolled flow. Our study confirms all of 
these observations; it also extends them at several different levels. First, we pay close at- 
tention to a net efficiency by computing the net power gained (positive efficiency) or lost 
(negative efficiency) in the presence of wall-actuation. Second, we conduct much more 
detailed study of the influence of traveling waves on velocity fluctuations; this is done 
by a thorough analysis of the influence of the wave speed, frequency, and amplitude on 
receptivity of full three dimensional fluctuations. Third, we confirm all of our theoretical 
predictions in Part 2 of this study, and highlight remaining research challenges. 

Our presentation is organized as follows: in § [2j we formulate the governing equations 
in the presence of traveling wave wall-actuation. The influence of control on the nominal 
bulk flux and the nominal net efficiency is also discussed in this section. A frequency 
representation of the NS equations linearized around base velocity induced by traveling 
waves is presented in § [3j We further discuss a notion of the ensemble average energy 
density of the statistical steady-state and describe an efficient method for determining 
this quantity in flows subject to small amplitude traveling waves. In § [4j we employ 
perturbation analysis to derive an explicit formula for energy amplification. This formula 
is used to identify the values of wave frequency and speed that reduce receptivity of the 
linearized NS equations; we show that the essential trends are captured by perturbation 
analysis up to a second order in traveling wave amplitude. We also discuss influence of 
amplitude on energy of velocity fluctuations and reveal physical mechanisms for energy 
amplification. A brief summary of the main results along with an overview of remaining 
research challenges is provided in § [5] 

2. Steady-state analysis 

2.1. Governing equations 
Consider a channel flow governed by the non-dimensional incompressible NS equations 

u t = -(u- V)u - VP + (l/R c )Au + F, = Vu, (2.1) 




with the Reynolds number defined in terms of the centerline velocity of the parabolic 
laminar profile U c and channel half- height 5, R c = U c S/u. The kinematic viscosity is 
denoted by the velocity vector is given by u, P is the pressure, F is the body force, V 
is the gradient, and A = V • V is the Laplacian. The spatial coordinates and time are 
represented by (x, y, z) and f, respectively. 

In addition to a constant pressure gradient, P X1 the flow is exposed to a zero- net- 
mass-flux surface blowing and suction in the form of a streamwise traveling wave (see 
figure [I] for illustration). In the absence of the nominal body force, F = 0, base velocity 
115 = (£/, V, W) represents the steady-state solution to \2.1\ subject to 



V{y = ±1) = T2acos(^(x - c*)), F = 0, 
U(±l) = Vy(±l) = W(±l) = 0, P- x = -2/R c , [ ' ] 

where u x , c, and a, respectively, identify frequency, speed, and amplitude of the traveling 
wave. Positive values of c define a DTW, whereas negative values of c define a UTW. 
The time dependence in V(±l) can be eliminated by the Galilean transformation, (x = 
x — ct, y = y, z = z, t = t). This change of coordinates does not influence the spatial 
differential operators, but it transforms the time derivative to di — dt — cd X: which adds 
an additional convective term to the NS equations 

u t = cu x - (u-V)u - VP + (l/i? c )Au + F, = V-u. (2.3) 

In new coordinates, i.e. in the frame of reference that travels with the wave, the wall- 
actuation ( |2.2| ) induces a two dimensional base velocity, 115 = (U(x, y), V(x, y), 0), which 
represents the steady-state solution to ( |2.3| ). Note that the spatially periodic wall actu- 
ation, V(y = ±1) = ^2 a cos (lj x x), induces base velocity which is periodic in x. 

The equations describing dynamics (up to a first order) of velocity fluctuations v = 
(ix, w) around base velocity, 115, are obtained by decomposing each field in ( |2.3[ ) into 
the sum of base and fluctuating parts, i.e., {u = 115 + v, P = P -hp, F = + d}, and by 
neglecting the quadratic term in v 

v £ = cv x - (ii6-V)v - (v.V)uft - Vp + (l/ii c )Av + d, = V-v. (2.4) 



Note that the boundary conditions (2.2) are satisfied by base velocity and, thus, velocity 



fluctuations acquire homogeneous Dirichlet boundary conditions. 

2.2. Base flow 

Let us first consider a surface blowing and suction of a small amplitude a. In this case, a 



weakly nonlinear analysis can be employed to solve (2.3) subject to (2.2) and determine 



the corrections to base parabolic profile; similar approach was previously used by |Jo- 



vanovic et al. 



(2006); Hcepffner & Fukagata (2009). Up to a second order in control 
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Figure 2. Second order correction to the nominal flux, Ub,2(c,uj x ), for (a) upstream waves; 
and (b) downstream waves in Poiseuille flow with R c = 2000. Note: the level sets are obtained 
using a sign-preserving logarithmic scale; e.g., 5 and —3 should be interpreted as Ub,2 = 10 5 
and Ub,2 = — 10 3 , respectively. 



amplitude a, U(x,y) and V(x,y) can be represented as 

U(x,y) = U (y) + aUtfay) + a 2 U 2 (x,y) + 0(a 3 ), 
V[x,y) = aV!(x,y) + a 2 V 2 (x,y) + 0(a 3 ), 

where Uo(y) = 1 — y 2 denotes base velocity in Poiseuille flow and (see Appendix [a| 

Ui(x,y) = Ui,- 1 (y)e- ,u - x + Ui, 1 (y)e iu ' x , 

Vx(x,y) = V 1 ,- 1 (y)e- iu ' x + V 1 , 1 (y)e iu ' x , 

U 2 (x,y) = U 2 ,o(y) + U 2 ,-2(y) e ~ 2i ^ x + U 2a {y)e 2 ^ x , 

V 2 (x,y) = y 2 ,_ 2 (y)e- 2i ^ + V 2 , 2 (y)e 2i "**. 



(2.5) 



HcepfTner & Fukagata (2009) recently showed that, in the absence of driving pressure 



gradient, the traveling waves induce nominal bulk flux (i.e., pumping) in the direction 
opposite to the direction in which the wave travels. While the first order of correction 
to the base velocity is purely oscillatory, the quadratic interactions in the NS equations 
introduce mean flow correction U2 J o(y) a ^ the- level of a 2 . The nominal bulk flux is deter- 
mined by Ub = (1/2) f\ U(y)dy where the overline denotes averaging over horizontal 
directions. In the presence of a pressure gradient, the nominal flux in flow with no control 
is Ub,o — (1/2) Uo(y) dy = 2/3, and the second order correction (in a) to Ub is given 

by Ub, 2 = (1/2) ^2,0(2/) dy. Figure 2 shows Ub, 2 as a function of wave frequency, oj x , 
and wave speed, c, in Poiseuille flow with R c = 2000. Except for a narrow region in the 
vicinity of c = 0, the upstream and downstream waves increase and reduce the nominal 
flux, respectively. Furthermore, for a given wave speed c, the magnitude of the induced 
flux increases as the wave frequency is decreased. 

Figure |3] is obtained by finding the steady-state solution of (2.3) subject to (2.2) using 
Newton's method. Originally, we have used base flow resulting from the weakly nonlinear 
analysis to initialize Newton iterations; robustness of our computations is confirmed using 
initialization with many different incompressible base flow conditions. The nominal flux 
and its associated nominal drag coefficient for a UTW with c = — 2 and uo x = 0.5, and 
a DTW with c = 5 and uj x = 2 are shown in this figure. The flux and drag coefficient 
of both laminar and turbulent flows with no control are also given for comparison. The 
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Figure 3. (a) The nominal flux, Ub(ol)\ and (b) the nominal skin-friction drag coefficient, 
Cf(ot), for a pair of UT Ws a nd a pair of DT Ws in Poiseuille flow with R c = 2000. The results 
are obtained by solving ( |2.3[ ) subject to (2.2), in the steady-state, using Newton's method; Ub 
and Cf of the uncontrolled laminar and turbulent flows are also shown for comparison. 



nominal skin-friction drag coefficient is defined as ( McComb||1991 ) 

C f = 2r w /U 2 B = -2P X /U 2 B , 

where r w is the nondimensional average wall-shear stress. For the fixed pressure gradient, 
P x = —2/R c , the nominal skin- friction drag coefficient is inversely proportional to square 
of the nominal flux and, in uncontrolled laminar flow with R c = 2000, we have Cf = 
4.5 x 10 -3 . The UTWs produce larger nominal flux (and, consequently, smaller nominal 
drag coefficient) compared to both laminar and turbulent uncontrolled flows. On the 
other hand, the DTWs yield smaller nominal flux (and, consequently, larger nominal 
drag coefficient) compared to uncontrolled laminar flow. In situations where flow with no 
control becomes turbulent, however, the DTWs with amplitudes smaller than a certain 
threshold value may have lower nominal drag coefficient than the uncontrolled turbulent 
flow; e.g., for a DTW with c = 5 and uo x = 2 this threshold value is given by a = 0.16 
(cf. figure [3(b)| . 

2.3. Nominal net efficiency 

For the fixed pressure gradient, the difference between the flux of the controlled and the 
uncontrolled flows results in production of a driving power (per unit horizontal area of 
the channel) 

Tlprod = —2P x (Ub,c — Ub,u), 

where Ub, c an d Ub, u are the nominal flux of the controlled and uncontrolled flows, 
respectively. On the other hand, the required control power exerted at the walls (per 
unit horizontal area of the channel) is given by ( |Currie||2003l ) 



VP\ 



VP 



(2.6) 



The control net efficiency is determined by the difference of the produced and required 
powers (Qua drio fc Ricco| |2004) 



n ne t — n 



prod 



EL 



where H ne t signifies the net power gained (positive II ne t) or lost (negative II ne t), in the 
presence of wall- actuation. 
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For small control amplitudes, the produced power can be represented as 

Uprod = Tlprod,0 + ^ ^prod,2 + 0(a 4 ), 

where 

The nominal required control power can be determined from (2.6) by evaluating the hor- 
izontal average of the product between base pressure, P, and base wall-normal velocity, 
V, at the walls. Since, at the walls, the nonzero component of V contains only first har- 
monic in x (cf. (2.5)), we need to determine the first harmonic (in x) of P to compute 
Hr Qq . Base pressure can be obtained by solving the two dimensional Poisson equation 

Pxx H - Pyy 

- (U X U X + 2V X Uy + Vy Vy) , (2.7) 

where P satisfies the following Neumann boundary conditions 

P V \ y=±1 = ((V XX + Vyy)/R C + CV X )\ y = ±1 . 

These are determined by evaluating the ^-momentum equation at the walls. For small 
values of a, weakly nonlinear analysis, in conjunction with the expressions for U and V 
given in § 2.2, can be employed to solve (2.7) for base pressure 

P(x,y) = aP^x.y) + 0(a 2 ), 
P^x.y) = P^-i^e-^ + P 1A (y)e [ "* x , 

where Pi,-i and P\ \ are determined from 

P",±i(y) ~ "iPiMv) = T2ioj x V li±1 (y)U^(y), 

Pf,_i(±l) = (^_i(±l) - ^V^^il))/^ + ci^-iCil), 
*m(±1) = (^m(±1) " oj 2 x V 1a (±1))/R c - ciou x V 1A (±l). 

Here, the prime denotes the partial derivative with respect to y, and the required power 
can be represented as 

U req = a 2 Il reqy2 + (9(a 4 ), 

iW = (Pi,-iV 1A + p m ^i,-i)U_i - (Pi,-iVi,i + Pi,i^i,-i)| y=1 . 

Since the second order correction to the nominal produced power, U pro d,2, is directly 
proportional to Ub,2, Tl P rod,2 is positive for UTWs and negative for DTWs. It turns out 
that smaller choices of uo x result in larger produced (for UTWs) or lost (for DTWs) power. 
One of the main points of this paper, however, is to show that it may be misleading to 
rely on the produced power as the only criterion for selection of control parameters; in 
what follows, we demonstrate that the required control power as well as the dynamics of 
velocity fluctuations need to be taken into account when designing the traveling waves. 

2.4. Nominal efficiency of laminar controlled flows 

We next examine the nominal efficiency of laminar controlled flows. Since we are inter- 
ested in expressing the nominal efficiency relative to the power required to drive flow with 
no control, we provide comparison with both laminar and turbulent uncontrolled flows. 
The net efficiency in fraction of the power required to drive the uncontrolled laminar flow 
is determined by 

%U net = U net /U = -a 2 \7r 2 (R c ;c,u x )\ + 0(a 4 ), (2.8) 

where IIo = — 2P X Ub,o and 7T2 = (H P rod y 2 — n reg) 2)/IIo. It can be shown that the second 
order correction to %H net , 7T2, is negative for all choices of c and uj x (see figure [4|. This 
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Figure 4. Second order correction to the nominal net efficiency, 7r2(c, cj x ), for (a) upstream 
waves; and (b) downstream waves in Poiseuille flow with R c = 2000. Note: the level sets are 
obtained using a sign-preserving logarithmic scale; e.g., —4 should be interpreted as 7T2 — — 10 4 . 
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Figure 5. The steady-state net efficiency, %n ne t, of laminar controlled flows as a function of 
control amplitude a for a UTW with (c = — 2, cj x — 0.5) and a DTW with (c = 5, u x = 2) 
at R c = 2000. The results are obtained by assuming that the uncontrolled flow (a) remains 
laminar; and (b) becomes turbulent. 



is because the required power for maintaining the traveling wave grows faster than the 
produced power as a is increased. In addition, figure [4] shows that \tt2\ is minimized for 
small wave speeds and for uj x G (1, 4). Formula (|2.8|) demonstrates that the control net 



efficiency is negative whenever the uncontrolled flow stays laminar (cf. figurej5(a)J). This 

"^([200 



is a special case of more general results by[Bewley|([2009j) and |Fukagata et a/.| ( |2009| ) which 
have established that any transpiration-based control strategy necessarily has negative 
net efficiency compared to the laminar uncontrolled flow. 

On the other hand, the net efficiency of the laminar controlled flow in fraction of the 
power required to drive the uncontrolled turbulent flow is determined by 



%n ne t — 



U B ,o 



( 



n 



turb 



U B 



turb 



u B . 



\ 



turb 



U B ,o 



a 



\7T 2 (RciC,bJ x 



+ 0(a 4 ), (2.9) 



>o 



where II 



turb 



-2P X U B 



turb' 



Since the bulk flux of the uncontrolled turbulent flow is 
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Figure 6. The wave amplitude, awx, for which the nominal net efficiency, %U ne t, is positive at 
R c = 2000 for (a) a pair of UTWs; and (b) a pair of DTWs. T he sol id curves are computed using 
Newton's method, and the dotted curves are computed using (2.10). The results are obtained by 



assuming that the controlled flow stays laminar while the uncontrolled flow becomes turbulent 



smaller than that of the uncontrolled laminar flow (i.e., Usjurb < Ub,o), it is possible to 



obtain a positive net efficiency for sufficiently small values of a. Note that formula (2.9) 
is derived under the assumption that the controlled flow stays laminar while the un- 
controlled flow becomes turbulent. Clearly, this formula represents an idealization since 
it assumes that laminar flow can be maintained by both UTWs and DTWs even with 
infinitesimal control amplitudes. It also indicates that increasing the control amplitude 
always decreases the nominal net efficiency. In a nutshell, the control amplitude needs to 
be large enough to maintain a laminar flow but increasing the control amplitude beyond 
certain value brings the efficiency down and eventually leads to negative efficiency. If the 
efficiency is negative, maintaining a laminar flow does not lead to any net benefit in the 
presence of control. This is further illustrated in figure |5(b)| where Newton's method is 
used to show that a positive net efficiency can be achieved for control amplitudes smaller 
than a certain threshold value (e.g., a < 0.05 for the DTW with c = 5 and uo x = 2). In 
addition, the net efficiency monotonically decreases as a is increased, as predicted by the 



weakly nonlinear analysis up to a second order in a (cf. (2.9)). 

An estimate for the maximum value of a for which a positive net efficiency is attainable 
can be obtained by solving the following equation (obtained using weakly nonlinear 
analysis) 

(1 - U B ,turb/U Bi o) ~ «Lx \ir 2 (Rc,c,uj x )\ = 0. (2.10) 
Figure [6] shows a max as a function of u x for different values of c. The dotted curves 



denote the approximation for a max obtained using (2.10). The values of a max (solid 
curves) obtained using Newton's method are also shown for comparison; we see that 
the predictions based on the second order correction capture the essential trends and 
provide good estimates for a max (especially for large wave speeds and for wave frequencies 
between 0.1 and 10). Figures [5] and [6] are obtained by assuming that the flow with 
control stays laminar while the flow with no control becomes turbulent. Whether or 
not the traveling waves can control the onset of turbulence depends on the velocity 
fluctuations; addressing this question requires analysis of the dynamics, which is a topic 
of § [3] and § |1J where we examine receptivity of velocity fluctuations around UTWs and 
DTWs to stochastic disturbances. 
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Figure 7. A Bloch wave d(x, y, z, t) defined in (3.2 ) is obtained by the superposition of weighted 
Fourier exponentials with frequencies (6 n ,k z ) | n ez, with weights determined by d n (y,k z ,t). 



3. Dynamics of fluctuations around traveling waves 

3.1. Evolution model with forcing 



A standard conversion of (2.4) to the wall- normal velocity (v) /vorticity (77) formulation 
removes the pressure from the equations and yields the following evolution model with 

forcing 

Eij) t {x,y,z,t) = Fi/>(x,y,z,t) + Gd(x,y,z,t), 

v(x,y,z,i) = Ci/>(x,y,z,t). 

This model is driven by the body force fluctuation vector d = (di, ^3), which can ac- 
count for flow disturbances. We refer the reader to a recent review article flSchmid|2007P 



and a monograph ( |Schmid fc Henningson|[2b01[ ) for a comprehensive discussion explain- 
ing why it is relevant to study influence of these excitations on velocity fluctuations. The 



internal state of (3.1 ) is determined by 1/? = (v, 77), with Cauchy (both Dirichlet and Neu- 



mann) boundary conditions on v and Dirichlet boundary conditions on 77. All operators 



in (3.1) are matrices of differential operators in three coordinate directions 7/, and z. 



Operator C in (3.1 ) captures a kinematic relation between i\) and v, operator G describes 
how forcing enters into the evolution model, whereas operators E and F determine in- 
ternal properties of the linearized NS equations (e.g., modal stability). While operators 
E, G, and C do not depend on base velocity, operator F is base-velocity-dependent and, 
hence, it determines changes in the dynamics owing to changes in U5 (see Appendix [B| . 
Moreover, for base velocity of § |2.2[ F inherits spatial periodicity in x from U5 and it can 
be represented as 

00 1 

F = F + J^ " 1 J2 e[rU!xXF l,n 

1 = 1 2 . 



where Fq and F^ r are spatially invariant operators in the streamwise and spanwise direc- 
tions and 1 signifies that r takes the values {—/,—/ + 2, — 2, /}. This expansion 
isolates spatially invariant and spatially periodic parts of operator F, which is well-suited 



for representation of (3.1) in the frequency domain. 



3.2. Frequency representation of the linearized model 

Owing to the structure of the linearized NS equations, the differential operators E, G, 
and C are invariant with respect to translations in horizontal directions. On the other 
hand, operator F is invariant in z and periodic in x. Thus, the Fourier transform in z can 
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be applied to algebraize the spanwise differential operators. In other words, the normal 
modes in z are the spanwise waves, e lkzZ , where k z denotes the spanwise wavenumber. 
On the other hand, the appropriate normal modes in x are given by the so-called Block 
waves ( Odeh fc Keller|1964| |Bensoussan et a/.|1978 ), which are determined by a product 
of e l0x and the 2tt/uj x periodic function in x, with G [0, lj x ). Based on the above, each 
signal in (|3.1|) (for example, d) can be expressed as 



d(x,y,z,i) 
d(x,y,k z ,t) 



e [k * z e [Gx d(x : y,k z ,t) 
d(x + 27r/u x ,y,k z ,t) 



k z G 



G [0, lj x ), 



where only real parts are to be used for representation of physical quantities. Expressing 
d(x,y,k z ,i) in Fourier series yields (see figure [7] for an illustration) 



d(x,y,z,t) 



d n {y,k z ,t) e 



\{0 n x + k z z) 



k z g R, G [0, u x ), 



(3.2) 



where {d n (y, k z ,t)} ne z are the coefficients in the Fourier series expansions of d(x, y, k z ,i). 

The frequency representation of the linearized NS equations is obtained by substituting 
(372b into JO) 



d t il> e (y,k z ,t) = Ae(k z )il> e (y,k z ,t) 
v e (y,k z ,t) = C e {k z )^ e (y,k z ,t). 



Be(k z )d (y,k z ,t), 



(3.3) 



This representation is parameterized by k z and and i^ e {y 1 k Zl t) denotes a bi- infinite 
column vector, tjj (y,k z ,t) = col{i/j(6 n ,y,k z ,t)} ne z. The same definition applies to 
de(y,k z ,i) and vo(y,k z ,t). On the other hand, for each k z and 0, Ae(k z ), Bo(k z ), and 
Ce{k z ) are bi-infinite matrices whose elements are one dimensional integro-differential 
operators in y. The structure of these operators depends on frequency representation of 



E, F, G, and C in (3.1). In short, Bo(k z ) and Co(k z ) are block-diagonal operators and 

00 

A e = Aoe + ^2 a 1 Aw, 



1 = 1 



where Aoq and Aw are structured operators (see Appendix |B| for more details). The 
particular structure of Aoe and Aw is exploited in perturbation analysis of the energy 
amplification for small control amplitudes a in § |3.4| 



3.3. Energy density of the linearized model 



Frequency representation (3.3) contains a large amount of information about linearized 



dynamics. For example, it can be used to assess stability properties of the base flow. 
However, since the early stages of transition in wall-bounded shear flows are not appro- 
priately described by the stability properties of the linearized equations (see, for example, 
Schmid & Henningson 2001 Schmid||2007|), we perform receptivity analysis of stochas- 



tically forced model (3.3) to assess the effectiveness of the proposed control strategy. 



Namely, we set the initial conditions in (3.3) to zero and study the responses of the lin- 



earized dynamics to uncertain body forces. When the body forces are absent, the response 
of stable flows decays asymptotically to zero. However, in the presence of stochastic body 
forces, the linearized NS equations are capable of maintaining high levels of the steady- 



state variance (Farrell & Ioannou 1993 Bamieh & Dahleh 2001 Jovanovic & Bamieh 



2005). Our analysis quantifies the effect of imposed streamwise traveling waves on the 
asymptotic levels of variance and describes how receptivity changes in the presence of 
control. We note that there are substantial differences between the problem considered 
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here and in |Jovanovic fc Bamieh ( 2005 ) ; these differences arise from lack of homogeneity 
in the streamwise direction which introduces significant computational challenges which 
we discus below. Furthermore, even though our study is similar in spirit to |JovanovIc" 
(2008), current work studies dynamics of fluctuations around spatially periodic base ve- 



locity, whereas | Jovanovic | (|2008 ) considered dynamics of fluctuations around time periodic 
base velocity. Theoretical framework for quantifying receptivity in these two conceptually 
different cases was developed by Fardad, Jovanovic & Bamieh (2008) and Jovanovic & 



Fardad (2008), respectively. 



Let us assume that a stable system (3.3) is subject to a zero-mean white stochastic 



process (in y and £), de(y,k z ,t). Then, for each k z and the ensemble average energy 
density of the statistical steady-state is determined by 

E(0,k z ) = lim (ve(',k z ,t),ve(',k z ,t)) 



trace ( lim £{vo(-,k z ,i) ®ve(',k z ,t)}) , 

\t — y CO / 



where (•, •) denotes the Li2\ — 1, 1] inner product and averaging in time, i.e., 



(V0, Vfl) 

£{v(.,t)} 



£ 



Vg(y,k z ,t) v e (y,k z ,t) dy }, 



lim — 



(3.4) 



v(-,t + r) dr, 



and v# (g) mq is the tensor product of v# with itself. We note that E(9,k z ) determines 
the asymptotic level of energy (i.e., variance) maintained by a stochastic forcing in (3.3). 
Typically, this quantity is computed by running DNS of the NS equations until the 



statistical steady-state is reached. However, for linearized system (3.3), the energy den- 
sity E(0,k z ) can be determined using the solution to the following operator Lyapunov 
equation QFardad et~aT|2008| ) 



A e (k z )Xe(k z ) + X e (k z )A* e (k z ) 



Be(k z )Bl(k z ), 



(3.5) 



E(0,fc*) = trace (M^) Q(k z ) C e (k z )) . 
Here, Xo(k z ) denotes the autocorrelation operator of i/) e , that is 

X e {k z ) = lim £ {il> e (-, k Z: t) (g) ^(-, k z , t)} . 

t — y co 

Since C|(fc 2 ) Ce(k z ) is an identity operator, we have 

CO 

E(6,k z ) = trace {X e {k z )) = J2 ^ce(X d (9 n ,k z j) , 



n — —co 



(3.6) 



where Xd(0 n ,k z ) denotes the elements on the main diagonal of operator Xq. We note 
that E also has an interesting deterministic interpretation; namely, if vo(-,k z ,t) denotes 
the impulse response of (3.3), then 



E(0, k z 



trace (v^(-, k z ,t) (g) v^(-, k z ,t)) dt. 



Thus, the same quantity can be used to assess receptivity of the linearized NS equations 
to exogenous disturbances of either stochastic or deterministic origin. 
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3.4. Perturbation analysis of energy density 



Solving (3.5) is computationally expensive; a discretization of the operators (in y) and 



truncation of the bi-infinite matrices convert (3.5) into a large-scale matrix Lyapunov 
equation. Our computations suggest that in order to obtain convergence of 

N 

E(0,k z ) * trace (X d (0 n ,fcs)), 

n = -N 

a choice of N between ten (for uo x ~ 0(1)) and a few thousands (for uo x ~ O(0. 01)) is 
required. Since we aim to conduct a detailed study of the influence of streamwise traveling 
waves on dynamics of velocity fluctuations, determining energy density for a broad range 
of traveling wave parameters, k z and still poses significant computational challenges. 
Instead, we employ an efficient perturbation analysis based approach introduced by |Far-| 



dad fc Bamieh| p008| ) for solving equation (3.5). For our problem, this approach turns 



out to be at least 20 times faster than the truncation approach. This method is well- 



suited for systems with small amplitude spatially periodic terms and it converts (3.5) 
into a set of conveniently coupled system of operator-valued Lyapunov and Sylvester 
equations. A finite dimensional approximation of these equations yields a set of algebraic 
matrix equations whose order is determined by the product between the number of fields 
in the evolution model (here 2, the wall- normal velocity and vorticity) and the size of 
discretization in y. While consideration of small wave amplitudes simplifies analysis by 
providing an explicit expression for energy density, it is also motivated by our earlier 
observation that large values of a introduce high cost of control which is not desirable 
from a physical point of view. 

It can be shown (see Appendix [C] for details) that the energy density of system (3.1) 
can be represented as 

oo 

E(0,k z ',R c ,a,c,u x ) = E o (0,k z ;R c ,uj x ) + ^a 21 E 2 i(6,k z ; R c ,c,u x ), < a < 1. (3.7) 

i = i 

Thus, only terms with even powers in a contribute to E, which in controlled flow depends 
on six parameters. Since our objective is to identify trends in energy density, we confine 
our attention to a perturbation analysis up to a second order in a. We briefly comment 



on the influence of higher order corrections in § 4.3 where it is shown that the essential 



trends are correctly predicted by the second order of correction. 



4. Energy amplification in Poiseuille flow with R c = 2000 

In this section, we study energy amplification of stochastically forced linearized NS 
equations in Poiseuille flow controlled with streamwise traveling waves. Equation (3.7) 
reveals the dependence of the energy density on traveling wave amplitude a, for < 
ft«l. However, since the operators in (3.3) depend on the spatial wavenumbers (0 and 
k z ), R C1 w x , and c, the energy density is also a function of these parameters. Finding the 
optimal triple (a,c,LJ x ) that maximally reduces the energy of the velocity fluctuations is 
outside the scope of the current study; instead, we identify the values of c and uj x that are 
capable of reducing receptivity in the presence of small amplitude streamwise traveling 
waves. Since we are interested in energy amplification of the transitional Poiseuille flow, 
we choose R c = 2000 in all of our subsequent computations. This value is selected 
because it is between the critical Reynolds number at which linear instability takes 
place, R c = 5772, and the value at which transition is observed in experiments and DNS, 
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R c w 1000. The same Reynolds number was used by Min et al. (2006) in their DNS 
study. 



4.1. Energy density of flow with no control 
We briefly comment on the energy density in uncontrolled Poiseuille flow with R c = 2000; 



for an in-depth treatment see Jovanovic & Bamieh ( 2005[ ). The appropriate normal modes 
in the uncontrolled flow are purely harmonic streamwise and spanwise waves, e lkxX e lkzZ , 
where k x denotes the streamwise wavenumber. Figure [8] illustrates the energy density of 
the uncontrolled flow as a function of k x and k z , which we denote by Eo(k x ,k z ). The 
streamwise constant fluctuations with 0(1) spanwise wavenumbers carry most energy in 
flow with no control. Namely, the largest value of Eo(k x ,k z ) occurs at (k x — 0, k z ~ 
1.78), which means that the most amplified flow structures (the streamwise streaks) are 
infinitely elongated in the streamwise direction and have the spanwise length scale of 
approximately 3.55, where S is the channel half-height. We note that these input-output 
resonances do not correspond to the least-stable modes of the linearized NS equations. 
Rather, they arise because of the coupling from the wall-normal velocity v to the wall- 
normal vorticity r\. Physically, this coupling is a product of the vortex tilting (lift-up) 
mechanism ( Landahi| [T975 ) ; the base shear is tilted in the wall- normal direction by the 
spanwise changes in i>, which lead to a nonmodal amplification of rj. This mechanism 
does not take place either when the base shear is zero (i.e., U' = 0), or when there are 
no spanwise variations in v (i.e., k z = 0). On the other hand, the least-stable modes 
(TS waves) of uncontrolled flow create a local peak in Eo(k x ,k z ) around (k z = 0, k x ~ 
1.2), with a magnitude significantly lower compared to the magnitude achieved by the 
streamwise constant flow stru cture s. Finally, we note that the uncontrolled energy density 
Eo(0,k z ;uj x ) as appeared in (3.7) can be obtained from Eo(k x ,k z ) using the following 
expression 

OO CO 

E (6,k z ]Lj x ) = E o (0 n ,M = ^ F>o(0 + nuj x ,k z ). 



In other words, for fixed uo x and 0, E (#, k z ;u x ) represents the energy density of velocity 
fluctuations that are composed of all wavenumbers k x = {6 + nu x } ne Z' In comparison, 
Eo(k x ,k z ) is the energy density of velocity fluctuations composed of a single wavenumber 
k x (see figure [9] for an illustration). 



Controlling the onset of turbulence by traveling waves: receptivity analysis 



15 




4.2. Energy amplification of flow with control 

We next consider energy amplification of velocity fluctuations in Poiseuille flow with 
R c = 2000 in the presence of both UTWs and DTWs. As shown in § |3]4j for small 
amplitude blowing and suction along the walls, the perturbation analysis yields an explicit 



formula for energy amplification (cf. (3.7)) 



^ ,a 7 v— = 1 + a g 2 {0,k z ]c,u x ) + ©(or 1 ), < a < 1. 

Thus, for small wave amplitudes the influence of control can be assessed by evaluating 
function g 2 = E2/E0 that quantifies energy amplification up to a second order in a. Sign of 
#2 determines whether energy density is increased or decreased in the presence of control; 
positive (negative) values of g 2 identify wave speed and frequency that increase (decrease) 
receptivity. Since function g 2 is sign-indefinite with vastly different magnitudes, it is 
advantageous to visualize g 2 using a sign-preserving logarithmic scale 

g 2 = sign(# 2 )logio(l + |^|). 

For example, g 2 = 5 or g 2 = —3, respectively, signify E2 = 10 5 Eo or E2 = — 10 3 Eo- 
Since g 2 (0, k z ; c,u x ) depends on four parameters, for visualization purposes, we confine 
our attention to cross-sections of g 2 by fixing two of the four parameters. We first study 
energy amplification of the modes with k z — 1.78 and k z — as a function of c and 
uj x ; these spanwise wavenumbers are selected in order to capture influence of control on 
streamwise streaks and TS waves, respectively. Since, in uncontrolled flow, streamwise 
streaks (respectively, TS waves) occur at k x = (respectively, k x = 1.2), fluctuations 
with 6 = (respectively, 0(u x ) = 1.2 — uo x \1.2/uo x \) are considered; these values of are 
chosen to make sure that streamwise streaks (respectively, TS waves) represent modes 
of the controlled flow as well. (Here, \a\ denotes the largest integer not greater than a.) 
We then analyze the energy amplification of disturbances with different values of and 
k z for a fixed set of control parameters c and uj x . Our analysis illustrates the ability of 
properly designed traveling waves to weaken the intensity of both most energetic and 
least stable modes of the uncontrolled flow. Direct numerical simulations of Part 2 show 
that this can be done with positive net efficiency. 
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Since most amplification in flow with no control occurs for fluctuations with (k x = 0, 
k z = 1.78), it is relevant to first study the influence of controls on these most energetic 
modes. In flow with control, the streamwise-constant flow structures are imbedded in the 
fundamental mode, i.e. fluctuations with = (cf. § 3.2). As the plots of ^(c, oo x ) m 
figures 10(a)| and 10(c)| reveal, the values of c and uj x determine whether these structures 
are amplified or attenuated by the traveling waves. Up to a second order in a, the 
control parameters associated with the blue regions in these two figures reduce the energy 
amplification of the uncontrolled flow. As evident from figure 10(a) , only a narrow range 
of UTWs with uj x < 0.1 is capable of reducing the energy amplification. However, since 
the required power for maintaining the nominal flow for such low frequency controls 
is prohibitively large (cf. figure |10(b)| ), the choice of UTWs for transition control is 
not favorable from efficiency point of view (receptivity reduction by these UTWs is 
further discussed in § 4.3). On the other hand, a large range of DTW parameters with 



Min et al. 



(2006) 



c > 1 and uj x > 0.1 is capable of making the controlled flow less sensitive to stochastic 
excitations (cf. figure |10(c)| ). Moreover, figure |10(d)| shows that the uj x > 0.1 region 
contains the smallest required power for sustaining the DTWs. These two features identify 
properly designed DTWs as suitable candidates for controlling the onset of turbulence 
with positive net efficiency (as confirmed by DNS in Part 2). 

It is noteworthy that traveling waves with parameters considered in 
(i.e., uj x = {0.5,1,1.5,2} and —4 < c < 0) increase amplification of the most energetic 
modes of the uncontrolled flow (cf. figure 10(a) ). This is in agreement with a recent study 
of |Lee et a/.| ( |2008 ) where a transient growth larger than that of the laminar uncontrolled 
flow was observed for UTWs with c = { — 1, — 2} and uo x = 1.5. Furthermore, it is shown 
in Part 2 that such UTWs promote turbulence even for initial conditions for which the 
uncontrolled flow stays laminar. 

The above analysis illustrates the ability of the DTWs to weaken the intensity of the 
most energetic modes of the uncontrolled flow; this is achieved by reducing receptivity to 
stochastic disturbances. However, an important aspect in the evaluation of any control 
strategy is to consider the influence of controls on all of the system's modes. In view of 
this, we next discuss how control affects the full three dimensional fluctuations. Since 
for a given uj x the energy amplification is symmetric around 6 = cj^/2, it suffices to 
only consider the modes with G [0,0^/2]. Figure 11 shows §2(0 \k z ) for a UTW with 
(c = —2, uo x — 0.5), and three DTWs with (c = 3, lj x = 1.5), (c = 5, uj x = 0.5), and 
(c = 5, uj x = 2). As evident from figure [l0| the selected UTW increases amplification of 
the fundamental mode with k z = 1.78; on the other hand, all three DTWs reduce energy 
amplification of modes with (0 = 0, k z = 1.78). Figure 11 further reveals that the largest 
change in amplification for all of these traveling waves takes place at (0 = 0, k z « 1.78), 
which is precisely where the uncontrolled flow contains most energy. This observation 
suggests presence of resonant interactions between the traveling waves and the most 
energetic modes of the uncontrolled flow. Additionally, as can be seen from figures ll(a)| 
and [TT(d)| the energy of modes with k z « is reduced by a UTW with (c = —2, uj x = 0.5) 
and a DTW with (c = 5, uj x = 2) for all 6. On the other hand, figure [11(b)] shows that a 
DTW with (c = 3, uj x = 1.5) increases amplification of fluctuations with (0.1 < 6 < 0.4, 
k z « 0); similarly, receptivity of fluctuations with (0.05 < < 0.45, k z « 0) is increased 
by a DTW with (c = 2, uj x = 0.5) (cf. figure 11 (c)[ ) . Thus, from the four considered 
cases, only a DTW with (c = 5, uo x = 0.5) can be used to inhibit intensity of full three 
dimensional velocity fluctuations (i.e., for all values of and k z ). 

While the fundamental mode is most influential in determining the effect of control 
on the energy amplification, figures [Tl(b) and 11(c) indicate that the modes with ^ 
and large spanwise wavelengths (i.e., k z ~ 0) can be significantly amplified by the 
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Figure 10. (a) and (c): Second order correction to the energy amplification, ^(c, uj x ), of the 
modes with (0, k z ) = (0, 1.78), in the presence of (a) UTWs; and (c) DTWs in Poiseuille flow with 
R c = 2000. (b) and (d): Second order correction to the nominal required power, Ii re g ; 2 ; c) , 
for (b) UTWs; and (d) DTWs. The dot and the square, respectively, denote (c = —2, uj x = 0.5) 
(as selected in |Min et a/.| ( |2006|) ) and (c = 5, uo x — 2). Note: the color plots are obtained using a 
sign-preserving logarithmic scale; e.g., #2 = 5 and §2 = —3 should be interpreted as E2 = 10 5 Eo 
and E2 = — 10 3 Eo, respectively. 



traveling waves. We thus take a closer look at how control affects the spanwise constant 
fluctuations. The TS waves are characterized by (k x = 1.2, k z = 0) and, for a given 
u x , they are imbedded in the modes of the controlled flow for fluctuations with 6(uo x ) — 
1.2 — uj x [1.2/uj x \. Figure 12 shows the second order correction g2(c,oj x ) to the energy 
amplification of the modes with k z = subject to both UTWs and DTWs. Note that 
figure [l2| correctly captures the increased intensity of the TS waves by DTWs with (c = 3, 
uj x = 1.5) and (c = 5, uj x = 0.5), as already observed in figures p"l(b)| and |l 1(c) We also 
see that the traveling waves considered in |Min etal. (2006) reduce energy of the TS waves 
(we recall that these promote amplification of the streamwise streaks; cf. figures |10(a)| 
and 12(a)). On the other hand, the DTW with c = 5 and uj x = 2 decreases energy 
amplification of both streamwise streaks and TS waves (cf. figures [I0(c)1 and [12(b)] ). The 
values of c and uo x capable of reducing the energy amplification (up to a second order in 
a) of both most energetic and least stable modes of the uncontrolled flow are marked by 
the dark region in figure Il3j 
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Figure 11. Second order correction to the energy amplification, g 2 (6,k z ), for traveling waves 
with (a) (c — —2, uj x — 0.5); (b) (c = 3, cj x = 1.5); (c) (c = 5, uj x — 0.5); and (d) (c = 5, cj x = 2) 
in Poiseuille flow with R c = 2000. 

4.3. Effect of control amplitude on energy amplification 

We next discuss influence of control amplitude on the energy amplification. We show that 
perturbation analysis (up to a second order in a) correctly predicts the essential trends. 
This is done by comparing perturbation analysis results with computations obtained 



using large-scale truncation of the operators in Lyapunov equation (3.5) 



The limit of the perturbation series (3.7) can be obtained by applying Shanks trans- 
formation ( Shanks||l955 Van Dyke| [l964) on the perturbation-analysis-based correction 
coefficients in ( |3.7| ). This transformation represents an effective means for providing con- 
vergence (respectively, faster convergence) to a divergent (respectively, slowly conver- 
gent) series (Sidi 2003). It turns out that Shanks transformation significantly increases 
the maximum value of a for which series (3.7) converges. Figure 14 shows the energy 



density of the fundamental mode, as a function of k z , in the uncontrolled Poiseuille flow 
with R c = 2000 and a flow subject to a DTW with (c = 5, uj x = 2, a = 0.025). The 



controlled flow results are obtained using truncation of series (3.7) up to a second order 



in a, and Shanks transformation up to a fourth order in a. Note that even though the 
second order correction overestimates the amount of receptivity reduction, it correctly 
captures the essential trends. 



Figure 15 compares energy density of the fundamental mode in uncontrolled Poiseuille 
flow with R c = 2000, and in the controlled flows subject to: (a) a UTW with c = — 2 
and uj x = 0.5, figure 15(a)| and (b) a DTW with c = 5 and uj x = 2, figure 15(b) The 
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Figure 12. Second order correction to the energy amplification, §2(0, lj x ), of the modes with 
(0(u x ), k z ) = (1.2 - u x [^J , 0), in the presence of (a) UTWs; and (b) DTWs in Poiseuille flow 
with R c = 2000. The dot and the square, respectively, denote (c = —2, uj x — 0.5) and (c = 5, 

LO x 2). 
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Figure 13. The dark regions identify values of wave speed and frequency that, up to a second 
order in a, suppress the energy amplification of both most energetic and least stable modes in 
Poiseuille flow with R c = 2000 subject to: (a) UTWs; and (b) DTWs. 



controlled flow results are obtained using Shanks transformation up to a fourth order in 
a, and they closely match the large-scale truncation results (hollow circles). Figure p~5(b)| 
shows that the properly designed DTWs with amplitudes equal to 5%, 10%, and 20% 
of the base centerline velocity reduce the largest energy density of the uncontrolled flow 
by approximately 28%, 60%, and 80%, respectively. It is noteworthy that substantial 
reduction is obtained at the expense of relatively small increase (compared to the laminar 
flow) in the nominal drag coefficient, which approximately increases by 1%, 4%, and 
13%. Further increase in the amplitude of a DTW with c = 5 and uo x = 2 results 
even in larger receptivity reduction. Part 2 demonstrates that this approach can be 
successfully used for controlling the onset of turbulence in flows subject to large initial 
disturbances. However, the power required for maintaining laminar flow under these 
conditions is prohibitively large, which limits the advantage of using DTWs for transition 
control from efficiency point of view. 

In contrast to DTWs, figure |15(a) demonstrates that the UTW with c = — 2 and 




Figure 14. Energy density, E(fe), of the fundamental mode — in Poiseuille flow with 
R c = 2000 and (c = 5, cjcc = 2, a = 0.025). The controlled flow results are obtained using 
perturbation analysis up to a second order in a, and Shanks transformation up to a fourth order 
in a. 




Figure 15. Energy density, E(k z ), of the fundamental mode = in Poiseuille flow with 
R c = 2000 subject to: (a) a UTW with c = -2 and u x = 0.5; and (b) a DTW with c = 5 
and uj x = 2. Shanks transformation up to a fourth order in a is used in computations. The 
truncation results (hollow circles) are obtained for (a) a = 0.016; and (b) a = 0.1. 




Figure 16. Energy density, E(k z ), of the fundamental mode = in Poiseuille flow with 
R c = 2000 subject to a UTW with c = —5 and uj x — 0.03. Shanks transformation up to a fourth 
order in a is used in computations. 
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uj x = 0.5 increases receptivity. We note that all of these trends are correctly captured 
by the second order correction (in a) to the energy amplification and that our results 



agree with the transient growth study of Lee et al. (2008). Furthermore, large energy 



amplification of the UTWs may be thought of as a precursor to flow instability; namely, 
it turns out that the UTWs destabilize the flow for a > 0.03 which is a smaller value 
compared to the amplitudes chosen in Min et al. ( 2006[ ) (a = 0.05 and a = 0.125, 
respectively). 
As described in 



4.2 



figure 10(a)| suggests that the UTWs with uj x < 0.1 can reduce 



the intensity of the most energetic modes of the uncontrolled flow. Here, we demonstrate 



that such UTWs lead to a very modest receptivity reduction. Figure 16 illustrates that 
a UTW with (c = —5, uo x = 0.03, a = 0.025) reduces energy amplification by about 
8%. On the other hand, modal stability analysis can be used to show that amplitudes 
as small as a « 0.03 make the flow linearly unstable. Therefore, relative to flow with no 
control, the UTWs at best exhibit similar receptivity to disturbances. 



For control amplitudes shown in figures [15] and [16j we have verified stability of fluctu- 
ations around base velocities in both UTWs and DTWs by computing the eigenvalues of 



the large-scale truncation of operator Ae(k z ) in (3.3). Compared to solving the truncated 



version of Lyapunov equation (3.5), perturbation analysis in conjunction with Shanks 
transformation provides much more efficient way for determining energy amplification. 
For example, while it takes four days on a PC to obtain the truncated results (hollow 
circles) in figure [l5(b)| the Shanks approximation is computed in four hours on the same 
PC. Moreover, once the correction coefficients in (3.7) have been determined, the energy 



amplification for a reasonably wide range of control amplitudes can be obtained at no 
further cost. 

4.4. Energy amplification mechanisms 

The energy of velocity fluctuations around a given base flow can also be obtained from 
the Reynolds-Orr equation (Schmid & Henningson | |2001| ). This equation can be used to 
elucidate the energy amplification mechanisms and facilitate better understanding of the 
influence of UTWs and DTWs on transitional channel flows. In this section, we consider 
the Reynolds-Orr equation for the fundamental modes (i.e., modes with 6 = 0; cf. equa- 
tion (3.2)). Our results reveal that, relative to the uncontrolled flow, the DTWs reduce 



the production of kinetic energy, thereby enabling the smaller receptivity to disturbances. 
As opposed to the DTWs, the UTWs increase the production of kinetic energy. For the 
streamwise-periodic base flow, 115 = (U(x,y), V(x,y), 0), the time evolution of the ki- 
netic energy of the fundamental modes, E(0 = 0, k z ;t) = {vo(-,k z ,t),vo(-,k z ,t)) \e = o, is 
governed by 



1 dE 

2 dt 



77 ^7 = -(u e ,U y ve) - (v e ,V y ve) - (ue,U x ue) - (v e ,V x u e ) + 

(l/Rc)((v e ,V xx v e ) + (v0,dyyV0) - k 2 z (V0,V0» + (V0,d0>, = 0, 



(4.1) 



where, for example, vq = co\{v(nuj x ,y,k z ,t)} ne z for the fundamental modes. In ( |4.1[ ), 
V xx is a diagonal operator with {(—nuj x ) 2 I} ne z on its main diagonal, / is the iden- 
tity operator, (•, •) denotes the L/2[—l, 1] inner product and averaging in time (cf. equa- 
tion (3.4)), and U y , V y , U x , and V x are block- Toeplitz operators whose rth sub-diagonals 
are determined by the rth harmonic in the Fourier series representation of U y (x,y), 
V y (x,y), U x (x,y), and V x (x,y) (see Appendix [B| for details). The first four terms on 



the right-hand-side of (4.1) denote the work of Reynolds stresses on the base shear and 
they contribute to production of the kinetic energy. The next group of terms represents 
viscous dissipation and the last term accounts for the direct work of the forcing on the 
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(d) uncontrolled, cj x = 0.5 (e) UTW, c = -2,uj x = 0.5 (f) UTW, a = 0.015 



Figure 17. Contribution of production and dissipation terms to energy density of = mode 
in Poiseuille flow with R c = 2000 subject to (a)-(c) a DTW with (c = 5, u x = 2); and (d)-(f) a 
UTW with (c = —2, uj x = 0.5). (a,d) uncontrolled flow; (b,e) second order corrections; and (c,f) 
controlled flows. In (c,f), the controlled flow results are obtained using approximation up to a 
second order in a, and the uncontrolled flow results are shown for comparison. 



velocity fluctuations. It can be shown that the direct work of d on v is balanced by a 
fixed portion of the viscous dissipation, and that the difference between the production 
terms and the remaining dissipation terms determines the energy density (cf. figure 17). 

In the steady-state limit, (4.1) can be used to obtain the following expression for the 
energy density of the fundamental mode, E(0, k z ) = lim^ _^oo E(0, k z \ t), 

E(0,fc*) = E p (0,k z ) + E d (0,fc*). 

Here, E p (0, k z ) denotes the contribution of production terms to the energy density and 
Ed(0, k z ) represents the joint contribution of viscous dissipation and the work of distur- 
bances 

E p (0,k z ) = -(R c /k 2 z ) ((u e ,U y ve) + (v , V y v e ) + (u e ,U x u e ) + (v , V x u e )) , 
E d (0,fc*) = (l/k 2 z )((v e ,V xx v e ) + (v G ,d yy v e )) + (R c /k 2 z ) (v,,d,), = 0. ' 2) 

In flows subject to small amplitude traveling waves, a perturbation analysis can be em- 
ployed to study the effect of each term on the right-hand-side of (4.2) on the energy 
density 

E(0,fc*) = E (0A) + a 2 E 2 (0,fc*) + (9(a 4 ), 
E (0A) = E , p (0,fcs) + E ,d(0,fc*), (4.3) 
E 2 (0,fc*) = E 2lP (0,fc z ) + E 2 , d (0,fc*), 
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where all the above terms can be readily determined from the solution to the Lyapunov 
equation (3.5). 

Figure 17 illustrates E p (0,k z ) and ~Ed(0,k z ) in the uncontrolled flow and in flows 
subject to a DTW with (c = 5, cj x = 2) and a UTW with (c = -2, uo x = 0.5). As 
expected, in the uncontrolled flow the joint contribution of dissipation and forcing is 
negative while the contribution of production is positive (see figures 
The energy density (solid curve) is determined by the sum of Eo, p and Eo,d, and it 
peaks at k z « 1.78. The second order corrections (in a) to E p and E^ are shown in 
figures [17(^)1 (for the DTW) and [l7(e)| (for the UTW). In flows subject to a DTW, the 
correction to E p is negative while the correction to E^ is positive. Furthermore, the effect 
of E2 jP dominates that of E2,d which implies that the DTW reduces the energy density 
of the uncontrolled flow (solid curve in figure 17(b) shows that a DTW introduces a 
negative correction to Eq). On the other hand, flows subject to a UTW exhibit opposite 
trends; the correction to E p is positive, the correction to E^ is negative, and since Ei2, p 
has the dominant eff ect, th e UTW increases the energy density of the uncontrolled flow 
(solid curve in figure 17(e) shows that a UTW introduces a positive correction to Eo). In 
figures 17(c) and 17(f)| perturbation analysis up to a second order in a is used to show 
E p and E d (solid curves) for a DTW with a = 0.025 and for a UTW with a = 0.015. 
Relative to the uncontrolled flow (symbols), the DTW decreases both production and 
dissipation terms. On the contrary, the UTW increases both of these terms. For both 
UTWs and DTWs, production dominates dissipation and determines whether the energy 
is increased or decreased. In addition, our computations show that {ueMy v e) is orders 
of magnitude larger than the other production terms. Moreover, (u, d yy u) completely 
dominates other dissipation terms. Therefore, the work of the Reynolds stress uv against 
the base shear U y is responsible for almost all of the energy production and the maximum 
viscous dissipation is associated with the wall-normal diffusion of the streamwise velocity 
fluctuation. These results are confirmed by DNS of the NS equations in Part 2. 



5. Concluding remarks 

This paper disentangles three distinct effects of blowing and suction along the walls on 
pumping action, required control power, and kinetic energy reduction. We have shown 
that analysis of dynamics is paramount to designing the streamwise traveling waves. 
If velocity fluctuations are well-behaved then the pumping action and required control 
power can be ascertained from the steady-state analysis. The proposed method uses 
receptivity analysis of the linearized NS equations to study the fluctuations' energy in 
transitional channel flows. Motivated by our observation that a positive net efficiency 
can be achieved by preventing transition, we develop a framework for design of the 
traveling waves that reduce receptivity to three dimensional body force fluctuations. 



Direct numerical simulations of the NS equations, conducted in Part 2 of this study (Lieu 
et a/.||20l0 ), verify that the traveling waves identified here are indeed an effective means 



for controlling the onset of turbulence. This demonstrates the predictive power of model- 
based approach to sensorless flow control; our simulation-free approach captures the 
essential trends in a computationally efficient manner and avoids the need for DNS and 
experiments in the early design stages. 

Our perturbation analysis has revealed that properly designed DTWs can significantly 
reduce energy amplification of three dimensional fluctuations, including the streamwise 
streaks and the TS waves, which makes them well-suited for preventing transition. The 
DNS of Part 2 confirm that transient response of fluctuations' kinetic energy can be main- 
tained at low levels using the values of wave frequency and speed that reduce receptivity 
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of the linearized NS equations. This facilitates maintenance of laminar flow; positive net 
efficiency can be achieved if the wave amplitude necessary for controlling the onset of 
turbulence is not prohibitively large ( Lieu et aL]|201Q ). On the other hand, we show that 
the UTWs are poor candidates for preventing transition for they, at best, exhibit sim- 
ilar receptivity to background disturbances as the uncontrolled flow. In particular, the 
UTWs considered by Min et al ( 2006 ) largely amplify the most energetic modes of the 
uncontrolled flow, thereby promoting turbulence even when the uncontrolled flow stays 
laminar ( Lieu et aZ7||201Q| ) . 

In spite of promoting turbulence, the UTWs may offer a viable strategy for reduc- 
ing skin-friction drag in fully developed turbulent flows. The DNS of Min et al (2006) 
and Lieu et al\ ( |2010[ ) suggest that the UTWs alter the dynamics of velocity fluctua- 
tions favorably in the turbulent regime, e.g., by reducing skin-friction drag coefficient 
compared to the uncontrolled flow. However, since the UTWs induce turbulent flow that 
departs from base flow obtained in the absence of velocity fluctuations, the model used 
in this work cannot be employed to explain utility of UTWs. This would require develop- 
ment of control-oriented models that contain essential physics of turbulent flows and, at 
the same time, are convenient for control design. For example, turbulent viscosity mod- 



els have been successfully used to determine the turbulent mean velocity (Reynolds 
Tiederman|1967 Reynolds fc Hussam||1972 ) and to identify the dominant turbulent flow 



structures (Del Alamo & Jimenez 2006; Cossu e^"al]|2009 Pujals et a/.|2009| ). Motivated 
by these successes, we intend to explore development of new control-oriented models that 
are capable of capturing the essential features of turbulent flow dynamics. 

The contribution of this paper goes beyond the problem of designing transpiration- 
induced streamwise traveling waves. The techniques presented here may also find use in 
designing periodic geometries and waveforms for maintaining the laminar flow or skin- 
friction drag reduction. Our work (i) suggests that strategies capable of reducing high 
flow sensitivity represent viable approach to controlling the onset of turbulence; and (ii) 
offers a computationally attractive (and simulation-free) method to determine the energy 
amplification of the linearized flow equations in the presence of periodic controls. 
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Appendix A. Base velocity 

In order to determine the corrections to base parabolic profile in flows subject to small 
amplitude traveling waves, we use a weakly nonlinear analysis to solve (2.3) subject 
to (2.2). We only present the equations for corrections up to a second order in a; similar 
equations can be obtained for higher order corrections. Stream functions, ^i } ±i(y), can be 
used to determine the first harmonic in Fourier series representation of the base velocity 
(cf. ( [2^1 )) 

Ui,±i(y) = *i,±i(j/), V li±1 (y) = TWx*i,±i(y), 
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where ^i y ±i(y) are solutions to 

(l/Rc)Al m * l9±1 ± iu x ((c - U )A LJx * 1 , ±1 + E#*i,±i) = 0, 
*i,_i(±l) = ±i/u x , tf M (±l) = TiMo *i,±i(±l) = 0. 

Here, A Ux = d yy —uo 2 with Dirichlet boundary conditions and = d yyy y — 2oj 2 dyy-\-oj^ 
with Cauchy boundary conditions. Moreover, C/2,0 is obtained by equating terms of order 
a 2 in the streamwise averaged x-momentum equation 

(l/RcM[ = F M - U^Vl^ + V l9 -!U[ tl - [/i,-!^, u 2 , (±i) = 0. 



Appendix B. Frequency representation of the evolution model 

We first describe how base velocity modified by the traveling waves enters in evolution 
model (3.1). Frequency representation of the evolution model is discussed next. It turns 
out that the components of base velocity determine coefficients of operator F in (3.1). 
For base velocity, u& = (U (x, y), V(x, y), 0), F is a 2 x 2 block-operator with components 



(1/R C )A 2 + ((AU) -{U- cI)A)d x - (AV)d y 
2V x d xy + U x (A-2d xx ) - (AV y ) + (2(AV)d x 
V X (A - 2d yy ) - 2U x d xy ) (d xx + d zz )~ x d xy , 
- (2(AV)d x + AV X + V X (A - 2d yy ) - 2U x d xy ) (a 

-(Uydt + VxidM + d,,)- 1 ^ 

(1/R C )A - (U x + (£/ 
where (fl^ + d zz )~ x is defined by 



AV X A 



F 12 
F 21 
F 22 



J yyz, 
d)d x -f 



J vv) 

Vdy) 



d zz )~ x d z 



V x {d xx ~\~ &zz) &xyi 



dzz)- L ■ f ^ 9 



Frequency representation (3.3) of the linearized evolution model (3.1) can be deter- 
mined using the following simple rules (Fa rdad et a/.|[2QQ8| ): 

(a) A spatially invariant operator L with Fourier symbol L(k x ) has a block-diagonal 
representation jCq = diag {L(0 n )} n e ^. For example, if L = c^, then = diag{i(# + 
nuj x )I} nG z- Operators G, C, To, and T/ ?r in (3.1) are spatially invariant and, thus, 
their representations are block-diagonal. 

(b) A spatially periodic function T(x) with Fourier series coefficients {T n } ne z has a 
^-independent block- Toeplitz representation 



T T_! T_ 2 
T x T T_! 
T 2 Ti T 



where the box denotes the element on the main diagonal of T. For example, T(x) = e _ir:c 
has a block- Toeplitz representation T = S r with the only non-zero element T_ r = I. 

(c) A representation of the sums and cascades of spatially periodic functions and 
spatially invariant operators is readily determined from these special cases. For example, 
a matrix representation of operator e~ irx d x is given by S r diag {i(0 + nuj x )I} n e %. 



T 



toep j- 



T, rri rri rri rri 
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Based on these, we get the following representations for Aq, Bq, and Co in (3.3) 



z = 1 



i = i 



B e = £ G 1 Ge, Qe = diag {G(<9 n )} nGZ , C e = diag {C(# n )} nGZ , 

where we have used the fact that Eq = diag {E(0 n )} n e % is an invertible operator. For 
convenience of later algebraic manipulations, we rewrite Aw as Aw = J2 r 2 / S- r Ai^e 
where A^ r e = diag {Ai jr (O n )} n e % = diag {E~ 1 (9 n+r ) Fi jr (O n )} neZ . In other words, for 
a given / > 1 operator Aw has non-zero blocks only on rth sub-diagonals with r G 
{— /, — / + 2, . . . , I — 2, /}. The frequency symbols of the operators E(Q n ), G(# n ), G(# n ), 
and F^ r (0 n ) are given by 

^o n (#nA) = (lARc)A 2 + i0 n (C# - (U - c)A), F 12 (# n A) = 0, 
F^(0 n ,k z ) = -ik z Ui F 22 (6 n ,k z ) = (1/R C )A - iO n (U ~ c), 



and 



l,r\ 



\6 n ((A ruJx U hr ) - Ui :r A - 2i(ru x )Vi, r dy) - 

((A roJx Vi t r) + VA)d y - iruj x U hr (A + 2<9 2 ) - wuj x (A ru)x U hr ) ■ 

(0 n /k 2 ) (26 n ((-A ruJx V hr )d y + iruj x Ui, r d yy ) - 

ruj x ((A roJx Vi yr ) + Vz >r (A - 2dyy))d y ) , 

(/c,/& 2 ) (2# n ((-A rW£c Vi, r ) + irc^, r ^) - 

rcj x ((A ru;x Vz, r ) + V hr (A - 2d yy ))) , 

-ik z (u[ r - (iruj x /k 2 ) V hr d y 



-Vl r dy 



rUi, r + i0 n ((\ruj x /k 2 )Vi,rd y - U^ r ) , 



where k 2 = # 2 + /c 2 A = - /c 2 and A r 



dyy — (juo x ) 2 with Dirichlet boundary 



conditions, and A = d yyyy — 2k d yy 



k with Cauchy boundary conditions. Operators 



E, G, and C are given by 
E(0 n ,k z ) = 



A 
/ 



, G(0 n ,k z ) = 



-\6 n d y 
ik z I 



-k 2 I 




-ik z d y 



-i<9 7 , 



C(9 ni k z ) 



KQn/k 2 )d y -ik z /k 2 
I 

i(k z /k 2 )dy W n /k 2 



Appendix C. Perturbation analysis of energy density 



As discussed in § 3.3 , steady-state energy density, E(0, fc^), of the linearized sys- 
tem (3.3), can be determined using the solution to the operator Lyapunov equation (3.5). 



For sufficiently small values of a, the solution of (3.5) can be expressed as a perturbation 



series Xq = ^m = o am ^mQ- After substituting into (3.5) and factoring out the terms 
with equal power in a, we have 



a 



a : Aoo Xoo + Xqq Aq 9 = —BqBq, 
Aoe X rn Q + X m Q Aq 9 = 



(CI) 

Z = l 
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Since operator Aoe is block-diagonal, X m e inherits the same structure as the right-hand- 
side of (CI). One can show that X m Q has non-zero blocks only on the first s < m odd 
(for odd m) or even (for even m) upper and lower sub-diagonals. Up to a second order 
in a, we have 

Xoo = Xo,oo, 

%2Q — 52^2,26* + <^2,06>, +^2,20^-2? 

where X m , s o = diag {X m ,s(6n)}n e z and S r is defined in Appendix |b| Substituting 
into (CI) yields 

Aoe X 0y0 o + <^o,o<9 Aq 9 = —BqBq, 
Aoe SiXi^io + SiXi^qAqq = — (Si Ai^-io X 0y oe + ^o,o<9 Al^S^j •> 

Aoe <^2,0<9 + <^2,0<9 Aqq = — ( ^4-2, 06> <^0,0<9 + <%,0<9 ^2, 0(9 + $-1 A\,\Q S\ X\,\q + 

S\ Ai-ie X^ 1Q S-\ + 5i A\_ ie <S_i + ^1*10 5-i ^1,10 <Si) • 

(C3) 

Finally, each block on the main diagonal of X m ^ s Q in (C3) is obtained from 



Ao(6 n )Xo,o{6n) + ^0,o(^n)^o(^n) 
A (fln-l)^l,l(fln) + X M (0 n )A5(0 n ) 
A (fln)^2,0(fln)+^2,0(e n )A5(fl n ) 



-B(0 n )B*(0 n ), 

— (Al : -i(0 n )X O: o(9 n ) + X 0? o(^n-l)^l,l(^n-l)) 5 

— (^2,o(#n)^0,o(#n) + -^0,0 (#n)^2,0 (^n) + 
^l,l(#n-l)^l,l(#n) + ^l,-l(#n+l)^i,l(#n+l) + 
^l,l(0n+l)^l,-l(#n+l) + ^1,1 (^n)^-l,l (On-l)) - 
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